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Abstract 

In this paper we consider an algorithmic technique more general than that 
proposed by Zharkov and Blinkov for the involutive analysis of polynomial ideals. 
It is based on a new concept of involutive monomial division which is defined for 
a monomial set. Such a division provides for each monomial the self-consistent 
separation of the whole set of variables into two disjoint subsets. They are called 
multiplicative and non- multiplicative. Given an admissible ordering, this sepa- 
ration is applied to polynomials in terms of their leading monomials. As special 
cases of the separation we consider those introduced by Janet, Thomas and Pom- 
maret for the purpose of algebraic analysis of partial differential equations. Given 
involutive division, we define an involutive reduction and an involutive normal 
form. Then we introduce, in terms of the latter, the concept of involutivity for 
polynomial systems. We prove that an involutive system is a special, generally 
redundant, form of a Grobner basis. An algorithm for construction of involu- 
tive bases is proposed. It is shown that involutive divisions satisfying certain 
conditions, for example, those of Janet and Thomas, provide an algorithmic con- 
struction of an involutive basis for any polynomial ideal. Some optimization in 
computation of involutive bases is also analyzed. In particular, we incorporate 
Buchberger's chain criterion to avoid unnecessary reductions. The implementa- 
tion for Pommaret division has been done in Reduce. 
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1 Introduction 



In modern times the Grobner bases method invented by Buchberger has become 
one of the most universal algorithmic tools for analyzing and solving polynomial equa- 
tions p, ||. Even in the general case, when the roots cannot be exactly computed, the 
method is still able to obtain valuable information about the solutions. In particular, 
it allows one to verify compatibility of the initial equations and compute the dimen- 
sion of the solution space. For the last few years notable progress has been achieved 
in extension of the Grobner bases method to non- commutative 0, Q and differential 
algebras 0, 0. 

On the other hand, already by the early 20s the foundation of a constructive ap- 
proach to algebraic analysis of partial differential equations was laid by Riquier and 
Janet ^ giving, among other things, answers to the same general questions of com- 
patibility and dimension. Later on, this approach, in the context of partial differential 
equations, was developed by Thomas and more recently by Pommaret |TI|. The 
main idea of the approach, as with the computation of a Grobner basis, is rewriting 
the initial differential system into another, so-called, involutive form . 

In the involutive approach, unlike the Grobner basis method, independent vari- 
ables for each equation are separated into two distinct groups called multiplicative 
and non-multiplicative. Such a separation is determined by the structure of the lead- 
ing derivative terms. A differential system is called involutive if its non-multiplicative 
derivatives are algebraic consequences of multiplicative ones. In doing so, Janet 
Thomas [1^ and Pommaret |Tl| used different separations of variables. 

Zharkov and Blinkov |T^, Q argued that the involutive technique along with the 
Grobner bases one can be used in commutative algebra. Based on Pommaret definition 
of multiplicative and non- multiplicative variables [ill], they proved, among other things, 
that an involutive basis is a Grobner one. Moreover, their computational experience 
demonstrated a reasonably high efficiency of the new algorithm when it terminates. The 
termination, however, does not hold, generally, for positive dimensional ideals, while 
for zero- dimensional ones it does for any degree-compatible monomial orderings ||14 |. 



Apart from that, the Pommaret involutive form of Grobner bases for zero-dimensional 
polynomial ideals reveals a number of rather attractive features 
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In the present paper we consider an algorithmic technique more general than that 
proposed in [0, [1^ for the involutive analysis of polynomial ideals. First of all, we 



introduce a new concept of involutive monomial division (Sect. 3) which leads to the 
self-consistent separation of the whole set of variables into multiplicative and non- 
multiplicative subsets. Given an admissible ordering, the separation is applied to poly- 
nomials in terms of their leading monomials. That concept generalizes the particular 
choice used by Janet H] , Thomas 
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and Pommaret for analysis of partial differen- 
tial equations. We characterize also important properties of noetherity, continuity and 
constructivity for involutive divisions (Sect. 4). Noetherity provides for the existence of 
a finite involutive basis for any polynomial ideal. The other two properties allows one 
to construct that basis algorithmically. It is shown that all the above three divisions are 
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continuous and constructive. Thomas and Janet divisions are also noetherian whereas 
Pommaret division is not. 

Given an involutive division, we define an involutive reduction and an involutive 
normal form (Sect. 5). As this takes place, we show that much like the Pommaret 
normal form, investigated in |T^, the general involutive normal form is also unique 
and linear. Then we define involutive systems by analogy with differential equations 
(Sect. 6). To be involutive, systems are required to satisfy the involutivity conditions, 
which form the basis for their algorithmic construction. 

We prove (Sect. 7) that any involutive basis, if it exists, is a special, generally 
extended, form of the reduced Grobner basis. Though it is unique for Pommaret 



division [|14[, generally, it may not be the it is shown by an explicit example. 

We propose an algorithm for construction of involutive polynomial bases (Sect. 8). Its 
correctness is proved for any continuous involutive division and for arbitrary admissible 
monomial ordering, while its termination holds, generally, for noetherian divisions. The 



algorithm is an improved and generalized version of one proposed in [Tj, [T^ , and has 



been implemented in Reduce for Pommaret division. The main improvement is the 
incorporation of Buchberger's chain criterion [|T6|. 



2 Preliminaries 

Let M = irf polynomial ring over the field K of characteristic zero. In 

this paper we use the notations: 

f,g,h,p,q are polynomials in R. 

a, b, c are elements in K. 

F, G, H are finite subsets of M. 

N is the set of non-negative integers. 

M = { xf^ ■ ■ ■ x^" I G N, 2 = 1, . . . , n } is the set of monomials in M. 

T = {au\ueM, aeK} is the set of terms in M. 

u, V, w, s, t are monomials or terms with nonzero coefficients. 

U, V, W are finite subsets of M. 

degi{u) is the degree of xi in u. 

deg{u) is the total degree of u. 

cf{f, u) E K is the coefficient of the term u of the polynomial /. 

Id{F) is the ideal in R generated by the polynomial set F. 

>- is an admissible monomial ordering with xi >~ X2 >~ ■ ■ ■ >~ Xn- 

lt{f) is the leading term of / w.r.t. the ordering >~. 

lc{f) = cf{f,lt{f)) is the leading coefficient of /. 

lm{f) = lt{f)/lc{f) is the leading monomial of /. 

lm{F) = { lm{f) I / G -F } is the set of the leading monomials of F. 

lcm{F) is the least common multiple of the set { lm{f) | / G -F }. 

If the monomial u divides the monomial v we shall write u\v. 



3 



3 Involutive Monomial Division 



Definition 3.1 We shall say that an involutive division L or L— division is given on M 
if for any finite set ?7 C M a relation \l is defined on t/ x M such that for any u,Ui & U 
and v,w & M the following holds: 

(i) . u\lw implies u\w. 

(ii) . u\lu for any u & U. 

(iii) . u\l{uv) and u\l{uw) if and only if u\l{uvw). 

(iv) . If u\lw and ui\lw, then u\lUi or ui\lu. 

(v) . If u\lUi and Ui\lw, then u\lW. 

(vi) . li V <^U and m G then mIlw w.r.t. U implies u\lw w.r.t. V. 

li u\l{w = uv), we say u is an involutive divisor of w, w is an involutive multiple of 
and V is multiplicative for m. In such an event we shall write w = u x v. If m is a 
conventional divisor of w but not an involutive one we shall write, as usual, w = u ■ v. 
Then v is said to be non-multiplicative for u. 

The conventional monomial division, obviously, satisfies condition (iv) only in the uni- 
variate case. The simplest bivariate example: x\{xy) and y\{xy) but -'x\y and -^y\x. 
Definition 3J- for each u eU provides separation of the set of variables 

{Xi, ...,Xn} = Ml{u, U) U NMl{u, U) 

into two disjoined subsets {Ml{u,U) fl NMl{u,U) = 0) of multiplicative Ml{u,U) 
and non-multiplicative NMl{u,U) variables. It is convenient to define an involutive 
division for a monomial set just by specifying the subsets of multiplicative and non- 
multiplicative variables to satisfy the conditions (iv)-(vi). The other conditions will be 
fulfilled by the construction. 

Given an involutive division L and a finite set U, for each u E U let L{u, f/) C M 
be the set of multiplicative monomials for u, that is, 

u\lV <^==^ vEuL{u,U). (1) 

Then it is easy to see that Definition admits another form: 



Definition 3.2 An involutive division L on M is given, if for any finite U G M and 
for any u eU there is given a submonoid L{u, U) of M satisfying the conditions: 

(a). If w G L{u, U) and v\w, then v G L(m, U). 
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(b) . liu,v eU and uL{u, U) n vL{v, U) ^ 0, then u E vL{v, U) oi v e uL{u, U). 

(c) . \iv eU and v e uL{u, U), then L{v, U) C L(u, U). 

(d) . If y C [/, then f/) C L{u, V) for all ueV. 



We consider three different examples of involutive division introduced by Janet ^ , 
Thomas ||T0| and Pommaret |]ll| for analysis of algebraic differential equations. In doing 
so, we give, firstly, the definition of multiplicative and non-multiplicative variables for 
each of the divisions, and, secondly, prove the fulfillment of the three extra conditions 



iv)-(vi) in Definition |3.1| equivalent to (b)-(d) in Definition ^ 



Definition 3.3 Thomas division ITOl. Given a finite set U , let 



hiiJJ) = max{ degiiu) \ u eU } . 

A variable Xi is considered as multiplicative ioi u E U if degi{u) 
multiplicative, otherwise. 



hi{U) and non- 



Definition 3.4 Janet division 0. Let f/ be a finite set. For each 1 <i <n divide U 
into groups labeled by non-negative integ ers dij . . . J di'. 

[di, . . . ,di] = { u G f/ I degj{u) = dj, I < j < i }■ 

A variable Xi is multiplicative for m G f/ if i = 1 and degi{u) = max{ degi{v) \ v E U }, 
or if i > 1, M G [di, . . . , di-i] and 

degi{u) = max{ degiiy) \ v G [di, . . . , }. 



Definition 3.5 Pommaret division [0. For a monomial xf^ ■ ■ ■x'^ with dk > ^ the 



variables Xj with j > k are considered as multiplicative and xj with j < k as non- 
multiplicative. For M = 1 all the variables are multiplicative. 

We note that 

• Thomas division does not depend on the ordering on the variables Xj. Janet 
and Pommaret divisions, as defined, are based on the ordering of the variables 
assumed in Sect. 2. 

• The separation of variables into multiplicative and non-multiplicative ones for 
Thomas and Janet divisions are defined in terms of the whole set U. Contrast- 
ingly, Pommaret division is determined in terms of the monomial itself, regardless 
of the others, and, by this reason, admits extension to infinite monomial sets, un- 
like Thomas and Janet divisions. 
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To distinguish the above divisions the related subscripts T, J, P will be used. 
Proposition 3.6 Thomas, Janet and Pommaret monomial divisions are involutive. 

Proof According to the above remark we must prove that the conditions (iv)-(vi) in 
Definition p.l| are satisfied. 

Let M G t/ be a Thomas divisor of w G M, that is, w = m x f . Then degi{v) = 
degi{w) — hi{U) if degi{w) > hi{U) and degi{v) = if degi{w) < hi{U). Thus, if w has 
an involutive divisor m, then w/u is uniquely defined, and, hence, u is unique in U. It 
implies also the property (v) for Thomas division, since u\tv for u, f G f/ if and only 
if u = V. The property (vi) also holds since any hi for V is less than or equal to the 
corresponding hi for U. 

Let now u,v & U he two different Janet divisors of w, such that degiiu) = degiiv) = 
di for 1 < i < k < n and assume, for definiteness, that degkiu) > deguiv). Then, 
since both u^v are members of the same group [c^i, . . . the variable Xk is non- 

multiplicative for V. Hence, if m is a Janet divisor of w such that degkiw) > degkiu) > 
degk{v), then v is not Janet divisor of w. In other words, similar to Thomas division, 
any monomial w G M cannot have different Janet divisors in any set U. A monomial 
group may only be decreased by diminishing the set U, which implies the relation (vi). 

Lastly, consider a Pommaret divisor u of the monomial w = Xi^ ■ ■ ■ x'^ with m < n. 
By definition, u constitutes a left subset of the string representation for w as it is shown. 

u 

W = 'Xi - -Xi - ■ ■■ ■ - Xm - ■ -Xm ■ (2) 

V ' V ' 

dx dm 

It makes evident the fulfillment of the conditions (iv) and (v) for Pommaret division 
while the condition (vi) trivially holds since the division does not depend on the set U 
at all. □ 

Proposition 3.7 For any finite set U and for any u & U , the inclusion Mt{u, U) C 
Mj{u,U) and, respectively, NMj{u,U) C NMt{u,U) holds. 

Proof If Xi G Mj{u, U), u E [di, . . . , then, by definition, 

degiiu) = max{ degiiv) \ v G [di, . . . , di^i] } < max{ degiiv) \ v E U } . 

Hence, Xi G Mxiu, U) implies Xi G Mjiu, U). □ 

Definition 3.8 A set U is called involutively autoreduced with respect to division L 
or L— autoreduced if it does not contain elements L— divisible by other elements in U. 

Proposition 3.9 If U is L— autoreduced, then any monomial w EM has at most one 
L— involutive divisor in U . 
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Proof This follows immediately from the property (iv) of involutive division. In terms 
of Definition |3]^ it means that uL{u, U) fl vL{y, f/) = for all distinct u^v E U /liU is 
involutively autoreduced. □ 



Proposition 3.10 / (7^/ . If a set U is autoreduced with respect to Pommaret division, 
then for any u e U Mp{u, U) C Mj(u, U) and NMj{u, U) C NMp{u, U), respectively. 

Proof Let u = xf^ ■ ■ ■ xf.'' G M be a monomial with dk > and v E U he its 
Pommaret divisor. Then, as follows from the representation (0), 
with 1 < m < k and 1 < r < dm- It means that v G [di, . . . ,dm-i]- Since U is 
autoreduced by Pommaret division, there are no other members of the same group 
with degree in Xm higher than r. Therefore, v is also a Janet divisor of u, and u/v, 
being Pommaret multiplicative for v, is also Janet multiplicative. □ 

Example 3.11 U = {xy,y'^,z} {x y y y z). 



monomial 


Thomas 


Janet 


Pommaret 


Mt 


NMt 


Mj 


NMj 


Mp 


NMp 


xy 


X 


y,z 


x,y,z 




y,z 


X 




y 


X, z 


y,z 


X 


y,z 


X 


z 


z 


x,y 


z 


x,y 


z 


x,y 



4 Involutive Monomial Sets 

Definition 4.1 Given an involutive division L, a set U is called involutive with respect 
to L or L— involutive, if any multiple of some element m G f/, is also (L—) involutively 
multiple of an element f G f/, generally, different from u. It means that 

(Vm G U) (Vw G M) (3f G U) [ v\l{uw) ] (3) 

or, in accordance with (|l|) and Definition |3.2| , 



Definition 4.2 We shall call the set VJu^uuM. the cone generated by U and denote it 
by C{U). The set Uu^^u u L{u,U) will be called the involutive cone of U with respect 
to L and denoted by Cl{U). 

Thus, the set U is L— involutive if and only if its cone C{U) coincides with its involutive 
cone Cl{U). 
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Definition 4.3 A finite L— involutive set t/ C M will be called L— completion of a set 
[/ C [/ if C{U) = C{U). If there exists an L— completion U of the set U, then the 
latter is said to be finitely generated with respect to L. An involutive division L is 
called noetherian if every finite set U is finitely generated. 



Proposition 4.4 Given a noetherian involutive division L, every monomial ideal U 
has a finite involutive basis. 



Proof This is an immediately consequence of Definition O and Dickson's lemma |0 . 
□ 

Proposition 4.5 Thomas and Janet divisions are noetherian. 

Proof Given a finite set U, consider the monomial h where, as given in 

the definition of Thomas division, hi = max{ degi{u) \ u E U }, and form the finite set 
C M of all the different monomials v such that v\h and u\v for some u E U. The set 
V, which contains, in particular, the monomial h and the initial set U, is involutive for 
Thomas division. Indeed, let w = ■ ■ ■ xf^ be a multiple of some u G V. If w G V, 
then, obviously, w G Cr(l^). Otherwise, let {ciij, . . . , dj^.} {k < n) he the nonempty set 
which contains all the exponents di [1 < i < n) in w such that di-^ > hi-^, . . . , di,, > h^,,. 
Then there exists v eV satisfying 

W = VXi^ ■ • • Xi^ 

Since degi^iv) = hi-^, . . . ,degi,,{v) = hi,,, v is a Thomas involutive divisor of w, and, 
hence, w G Ct{V). 



Furthermore, from Proposition 3.7 it follows that there is a set of Vi C y which is 



a Janet completion of U. □ 

Definition 4.6 Multiplication of a monomial m G t/ by a variable x is called a pro- 
longation of u. Given an involutive division specified by the set U, the prolongation is 
called multiplicative if x is multiplicative for u and non-multiplicative, otherwise. 

In the construction of involutive sets the following concept of local involutivity plays 
the crucial role and admits the direct extension to polynomial sets (see Sect. 6). 

Definition 4.7 A set U is called locally involutive with respect to the involutive di- 
vision L if any non-multiplicative prolongation of any element in U has an involutive 
divisor in U, that is, 

(Vm G U) (Vx, G NMl{u, U)) {^v eU) [v\l{u- X,) ] (4) 
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In accordance with Definition [4.1| , the conditions (|^), apparently, are necessary for in- 
volutivity of U. Generally, however, they are not sufficient, as the next simple example 
shows. 



Example 4.8 Let L be an involutive division on M C K[x, y, z] defined by the table 



monomial 


M 


NM 


1 


x,y,z 




X 


X, z 


y 


y 


x,y 


z 


z 


y,z 


X 


ueM\ deg{u) > 2 




x,y,z 



It is easy to see that all properties listed in Definition ^TT] ( p.2| ) are satisfied, and the 
set U = {x,y, z} is locally involutive. For instance, x ■ y = y x x. However, U is not 
involutive since none -u G M with degxiu) > 0,degy{u) > 0,degzU > 0, e.g. xyz, has 
involutive divisors in U. 

The following definition and theorem enable one to reveal involutive divisions providing 
involutivity of every locally involutive set. 

Definition 4.9 An involutive division L will be called continuous if for any finite set 
U and for any finite sequence {ui}(i<i<fc) of elements in U such that 

(Vi < k) {3xj G NMl{u^, U)) [ ■ ] (5) 

the inequality Ui ^ uj for i ^ j holds. 

Theorem 4.10 // an involutive division L is continuous then local involutivity of any 
set U implies its involutivity. 

Proof Let set U be locally involutive, and such that any sequence in U satisfying (^) 
has no coinciding elements. We must prove that U satisfies (|^). Take any u E U 
and any G M and show that there is f G f/ such that v\l{uw). If u\l{uw) we are 
done. Otherwise, there is Xk^ G NMl{u,U) such that w contains Xk^- Then u ■ Xk^ 
has an involutive divisor vi G U. If Vi\l{uw) we are done. Otherwise, there are 
Xk2 G NMl{vi, U) and f 2 G f/ such that uw jvi contains and V2\L{yi • x^J. Going 
on, we obtain the sequence m, f 1, f2, . . . of elements in U satisfying (^). By construction, 
each element of the sequence divides uw. Since all the elements are distinct and uw 
has a finite number of distinct divisors, it follows that the above sequence in U is finite, 
and, hence, it ends up with an involutive divisor of uw. □ 

Corollary 4.11 Thomas, Janet and Pommaret divisions are continuous. 
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Proof Let U he a. finite set, and {Mi}(i<i<fc) be a sequence of elements in U satisfying 
the conditions @. We sliall sliow tliat tliere cannot be coinciding elements in the 
sequence for three divisions. 

It is ease to see that Ui+i\T{ui ■ XkJ imphes Uj+i = Ui ■ x^.. Indeed, suppose that 
Ui ■ Xk^ = Ui+i X Vi+i what means -iXk^\vi+i. If Vi+i would contain any other variable 
Xj^, then it would mean that degxj.{u.i) > degx^.^Ui+i), and, hence, Xj^ could not be 
multiplicative for Uj+i. Therefore, any Thomas sequence satisfying consists of 
distinct elements. 

If Ui+i\j{ui ■ XfcJ, then from definition of Janet division it follows that -Uj+i >~Lex Ui, 
where >- Lex is the lexicographical ordering corresponding to the choice of variable order 
xi >~ X2 >~ ■ ■ ■ >~ Xn a-^ assumed in Sect. 2. It is now obvious that Ui ^ uj for i j for 
Janet division. 

Let now Mj+i|p('Uj-Xfc.). Then the representation (||) shows clearly that Wj+i yRevLex 
Ui where y-RevLex is the reverse lexicographical ordering on M induced by the assumed 
variable order. □ 

With an eye to the below described algorithms based on examination of non-multiplica- 
tive prolongations only, we impose, in addition to continuity, one more requirement on 
an involutive division. 

Definition 4.12 We shall say that a continuous involutive division L is constructive 
if for any U C M, u E U, Xi E NMl{u, U) such that u ■ Xi ^ Cl{U) and 

(Vf E U) (yxj E NMiiv, U)) {v ■ Xj\u ■ Xi, V ■ Xj ^ u ■ Xi) [v- Xj E Cl{U)) ] (6) 

the following condition holds: 

iywECLiU)) [u-Xi^wL{w,UU{w})]. (7) 



Proposition 4.13 Thomas, Janet and Pommaret divisions are constructive. 

Proof Let T be Thomas division. Suppose there is ui E U, and v E T{ui, U) such 
that u ■ Xi = UiV X w, w E T{uiv, U U {uiv}). From Definition |3.3| it follows that if 
there exists Xj\w and ~^Xj\v for some 1 < j < Ti, then Xj E Mt{ui,U). This implies 
w E T{ui, U) and u ■ Xi E UiT{ui, U). 

Consider now Janet division J, and let m ■ Xj be a non-multiplicative prolongation 
which has no Janet divisors in f/, and for which the condition (^) holds. Assume for a 
contradiction that there \s Ui eU and v E J{ui, U) satisfying 

u ■ Xi = uiv X wi, Wi E J{uiV, U U {uiv}). 

Because v ^ 1 and wi ^ 1, select minimal j,m such that Xj\v and Xm\wi. It is 
easy to see that i < min{j,m}. Otherwise, by Definition |3.4|, we would have either 
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Xj ^ J{ui, U) if i > j or Xm ^ J{uiv, U U {uiv}) if z > m. Note that the equahty i = m 
impossible since U U {uiv} as well as any other monomial set is Janet autoreduced. 
Thus, Ml >~Lex u where >~Lex is the lexicographical ordering induced by the variable 
order Xi >~ ■ ■ ■ >- Xn- If monomial uxi is obtained by non-multiplicative prolongations 
of several elements in f/, then we suppose that u is lexicographically maximal from 
all of them. Since wi is non-multiplicative for mi, there is x^Jwi such that Xki € 
NMj{ui, U U {uiv}). Then, by condition (^), we can rewrite 

U-Xi = [Ui- XfcJ = [U2 X W2) = [U2 ■ XfcJ = ■ ■ ■ , 

where u -<Lex ^1 -<Lex U2 -<Lex ' ' ■• Continuity of Janet division implies termination 
of this chain with some ui & U such that u ■ Xi & UiJ{ui, U) what contradicts our 
assumption u ■ Xi ^ Cj{U). 

For Pommaret division condition follows directly from the property (v) in Def- 



inition 3.1 



□ 



Theorem 4.14 Let U he a non-involutive finitely generated set with respect to a con- 
structive division L. Then there is a procedure of completing U to an L—involutive set 
U D U based on enlargement of U by non-multiplicative prolongations of its elements. 



Proof Given U, by Definition |4.3| , there exists a finite L— completion U of U. We 
claim that U contains some non-multiplicative prolongations of elements in U. Assume 
for a contradiction that there are no such elements in U. Since set U is not involutive, 
there exist non-multiplicative prolongations of elements in U which have no L— divisors 
in U. 

Take any admissible ordering -< and select u eU with a non-multiplicative prolon- 
gation u ■ Xi which is not L— multiple of any element in U, and which is the lowest with 
respect to -<. Because U is involutive, and, by the above assumption, u ■ Xi ^ U, there 
is V E U \ U and 1 ~< w E M such that u ■ Xi = v x w, w E L{v, U). From the condition 
C{U) = Cl{U) it follows that v is multiple of some ui E U with degiui) < degiy). 

Show that V E Cl{U). If ui L— divides v, then we are done. Otherwise, there exists 
Xki\{v/ui), Xki E NMl{ui,U), and we can rewrite 

V = Ui = (Ml ■ XfcJ = {U2 X W2) = {U2 ■ Xk, 



until, by continuity of L, we come to an involutive divisor Um E U of v at some step of 
this rewriting procedure. This contradicts the constructivity condition (^, and, hence 
u ■ Xi E U. 

Now instead of U take Ui = U U {u ■ Xi} where Ui E U and Ui ■ Xi^ E U \ U . If 
set Ui is not involutive, then it can be further completed by the corresponding lowest 
non- multiplicative prolongation in Ui. Since the set U is finite, by repeating this 
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completion procedure, in a finite number of steps we construct tfie set f/ C f/ whicli is 
an L— completion of U. □ 

As an immediate consequence of the above described procedure of completing a set U 
by non-multiplicative prolongations of its elements we have the following corollary. 

Corollary 4.15 IfU is a finitely generated set with respect to a constructive involutive 
division, then there is the unique minimal involutive completion U of U such that for 
any other completion U the inclusion If U holds. 

The following algorithm, given a constructive division L, computes the minimal invo- 
lutive completion If for any finitely generated set U and any fixed admissible ordering 



-<. Its correctness and termination are provided by Theorem [4.14 

Algorithm InvolutiveCompletion: 

Input: U, a finite monomial set 
Output: U, an involutive completion of U 
begin 

U:=U 

while exist u eU and x G NMl{u, U) such that 
u ■ X has no involutive divisors in U do 
choose such u and x with the lowest u ■ x w.r.t. -< 
il ■.= UU{u-x} 

end 
end 



Example 4.16 (Continuation of Example p.ll|) . The minimal involutive bases of the 



set U = {xy,y'^,z) {x >~ y >~ z) for Thomas, Janet and Pommaret divisions are 

Ut = {xy, ^/^ z, xz, yz, xy^ xyz, y'^z, xy'^z} , 
Uj = {xy,y'^,z,xz,yz} , 

TT r 2 2 2 k 777 T 

Up = {xy, y ,z,xz,yz,x y,x z,...,x y,...,x z,...}, 
where /c, m G N. These bases can be easily derived from U using algorithm Involutive- 



Completion. Note that Uj C Ut and Uj C Up in agreement with Propositions 
and |3.10| . This example explicitly shows that Pommaret division is not noetherian. 
However, for another ordering z >~ y >~ x the set U is finitely generated, and then 
Up = U. 

5 Polynomial Reduction 



In this section we generalize the results obtained in |]T3|, |TJ] for Pommaret division to 
arbitrary involutive division. 
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Definition 5.1 Given a finite polynomial set F C M and an admissible ordering >-, 
the concept of multiplicative and non-multiplicative variables for / G -F is to be defined 
in terms of lm{f) and the leading monomial set lm{F). 

Therefore, as soon as we have polynomials rather than monomials, any involutive 
division is to be determined on the basis of some admissible ordering, even when it 
does not depend on the latter for the pure monomial with Thomas division. 

The concepts of involutive polynomial reduction and involutive normal form are 
introduced similar to their conventional analogues (Buchberger, 1985) with the use of 
involutive division instead of the conventional one. 

Definition 5.2 Let L be an involutive division L on M, and let F be a finite set of 
polynomials. Then we shall say: 

(i) . p is L— reducible modulo f ^ F ii p has a term t = au&T{a^O) such 

that u = lm{f) X t>, t> G L{lm{f) , lm{F)) . It yields the L— reduction p g = 
p- {a/lc{f)) f X V. 

(ii) . p is L— reducible modulo F if there exists f F such that p is L— reducible 

modulo /. 

(iii) . p is in L— normal form modulo F if p is not L— reducible modulo F. 

We denote an L— normal form of p modulo F by NFl{p, F). In contrast, a conventional 
normal form will be denoted by NF{p,F). As an involutive normal form algorithm 
one can use, for example, the following: 

Algorithm InvolutiveNormalForm: 

Input: p, F 

Output: h = NFl{p,F) 
begin 

h := p 

while exist f & F and a term u of h such that 

lm{f)\L{u/cf{h,u)) do 

choose the first such / 

h:=h-{u/ltif))f 
end 
end 

Correctness and termination of this algorithm can be proved, apparently, as they do 
for the conventional normal form algorithm |^. Since involutive reductions form a 
fixed subset of the conventional ones, generally, NFl{p,F) ^ NF{p,F). 

Definition 5.3 A set F is called involutively autoreduced with respect to the given 
involutive division L, or L— autoreduced, if the set lm{F) is L— autoreduced and every 
/ G F has no terms t = cf{f, t)u^ lt{f) with c/(/, t) ^ and m G CL{lm{F)). 
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Given an involutive division L and a finite set F, the following algorithm returns an 
L— autoreduced set H, denoted by if = Autor educe l{F) , and such that Id{F) = 
Id{H). 

Correctness of the algorithm is obvious from the while-loop structure. Since the 
underlying set of involutive interreductions is a subset of the conventional interreduc- 
tions, its termination follows from that for the conventional autoreduction 0, |^. 

Algorithm InvolutiveAutoreduction: 
Input: F 

Output: H = Autor educe l{F) 
begin 
H := F 

while exist h E H and g E H \ {h} 

such that h is reducible modulo g do 

choose the first such h 

H' ■.= H\{h} 

h' := NFL{h,H) 

if h' = Q then H := H' 

else H:=H'U{h'} 
end 
end 



Theorem 5.4 // set F is L— autor educed, then NFl{p,F) = if and only if p is 
presented in terms of a finite sum of the form 

peSpCR, 8F = {Y.fiXUi,\fieF,UijET} (8) 
with lm{uij) 7^ lm{uik) for j ^ k. 



Proof If NFl{p,F) = 0, then, by Definition ^]2| of involutive reductions, at 

each intermediate reduction step the current value p' of p is rewritten as p' — p" = 
p' — fi X Uij. Since the reduction chain is finite by admissibility of an ordering >~, the 
representation (||) holds. 

<^=: Let p is given by expression (^. Firstly, we show that lm{p) has an involutive 
divisor in the set lm{F). For this purpose select the leading term in the right hand side 
of (H). It has the form s = lt{fi x Uij) = lt{fi) x Uij with some i,j and cannot appear 
in any other term lt{fk) x um- Otherwise, the underlying monomial s/lc{s) would have 
two involutive divisors lm{fi) and lm{fk) what, by Proposition ^]9|, would contradict 
the involutive autoreduction of F. Secondly, since p is involutively reducible, after 
each reduction step the representation (||), obviously, still holds providing the further 
reductions until the chain stops when we obtain zero at a certain step. It just means 
that A^Fi(p,F) = 0. □ 
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Corollary 5.5 If set F is L—autoreduced, then the L— normal form, for an arbitrary 
algorithm of its computation and for any polynomials pi,P2 and p, has the properties: 

(i) . Uniqueness: if hi = NFl{p,F) and h2 = NFl{p,F) then hi = h2- 

(ii) . Linearity: NFl{pi +P2,F)= NFl{pi, F) + NFl{p2, F) . 



Proof (i) By an involutive normal form algorithm, hi = p — Y^^j fi x Uij and /i2 



P-Hij fi X Vi. 



Therefore, hi — h2 has the representation (|^), and NF^^hi — /12, -F) = 
On the other hand, since hi and /i2 are normal forms, they have no 
/i2. Hence, we have hi = h2- 



by Theorem pA 

involutive divisors and so does hi — 
(ii) Denote pi + p2 by and let 



hi = NFl{pi, F), h2 = NFl{p2, F), hs = NFLip3, F) . 



Then NFiiJi^ — hi — h2,F) = h^ — hi — /12, since none of hi,h2,h3 has involutive 



divisors in lm{F). In addition, because hk = Pk — fi ^ Vk^ij {k = 
hs — hi — h2 G Si?. Thus, by Theorem |^ , NF^^h^ — hi — h2,F) 
h3 = hi + h2. 



1,2,3), 



we have 



0, and, hence. 



□ 



6 Involutivity Conditions 

Definition 6.1 Multiplication of a polynomial / G -F by a variable x is called the 
prolongation of /. Given an involutive division specified by the set lm{F), the prolon- 
gation is called multiplicative if x is multiplicative for lm{f), and non-multiplicative, 
otherwise. 

Definition 6.2 An L—autoreduced set F is called (L—) involutive basis of Id{F) if 

(V/ G F) (Vn G M) [ NFlUu, F) = ] . (9) 



Proposition 6.3 Let F he an involutive polynomial basis. 
lm{F) is also involutive. 



Then the monomial set 



Proof It follows immediately from Definitions [4.1|, and 



□ 



It is clear from Definition |6.2| and the linearity of the involutive normal form, by Corol- 
lary ^.5| , that an involutive basis provides decision of the ideal membership problem. 
Hence, we have the following corollary. 
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Corollary 6.4 If set F is L—involutive, then p G Id{F) if and only if NFl{p, F) = 0. 
In this case, obviously, the equality Sp = Id{F) holds. 

The definition of involutive polynomial sets is the direct extension of that for involutive 
monomial sets in Sect. 4. The theorem below imparts the constructive characterization 
of involutivity, which is the heart of the involutive algorithms. 

Theorem 6.5 An L—autoreduced set F is involutive with respect to a continuous in- 
volutive division L if and only if the following conditions of local involutivity hold 

(V/ G F) (Vx, G NML{lm{f), lm{F))) [ iVFz.(/ • x„ F) = ] . (10) 



Proof Since Xi G M we are done. 

'^=: An immediate consequence of (p!0| ) is local involutivity of the set lm{F) in 



accordance with Definition f4.7| . Then, by continuity of division L, this set is involutive. 
Thus, for any f & F and any -u G M the monomial lm{f) ■ u has the involutive divisor 
Imiyg), g e F. 

We claim that the polynomial / ■ u can be presented as follows 

f . u = g X V + J2fi^ij ^ (11) 
ij 

where v, Vij G T, /j G -F and relation /m(/ ■ u) = lm{g x v) >- ImlfiVij) holds for any 
term of the sum. Indeed, if u is multiplicative for / we are trivially done. Otherwise 



u contains Xk G N Mi^f ,lm[F)). Then, the local involutivity of F, by Theorem |5.4| , 
yields the representation 

f ■ Xk = giy- ui + ^fiX Uij (12) 

ij 

with G -F and lm{f ■ Xk) = lm{giUi) >- fiUij for any term under the summation sign. 
If monomial u/xk is multiplicative for gi, then (|TT]) immediately follows from (T^) with 



g = gi and v = Uiu/xk- Otherwise, multiply both sides of (0) by u/xk, take a variable 
Xm € NML{gi,lm{F)), which is contained in u/xk, and apply the local involutivity 
conditions for gi ■ Xm- It gives the relation 

/ . M = (5(2 X U2)UiU/ {XkXm) + 51 fi^ij (^3) 

where inequality lm{g2)uuiU2/ (xkXm) >- l^ifiUij) holds for all i,j. If uui/{xkXm) is 
still non-multiplicative for g2 the relation ([l^) can be further rewritten by using the local 
involutivity conditions until we obtain relation ([TT]). This is guaranteed by continuity 
of involutive division L, because all the polynomials gi,g2, . . . G F are distinct, since 
their leading monomials, by construction, form the sequence satisfying (B|). 
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Next, similar rewriting the every term fiVij in ([TT]) gives fiVij = fkX Wk + T.im fi^im 
with lm{fiVij) = lm{fk x Wk) >- lm{fiWim)- Proceeding with this way, by admissibihty 
of ordering -<, we find, in a finite number of steps, that f ■ u & Sp. □ 

The next definition of partial involutivity is useful for the algorithmic construction of 
involutive bases as we show below. 

Definition 6.6 Given v G M and an L— autoreduced set F, if there exist f E F such 
that lm{f) -< V and 

(V/ G F) (Vm G M) (Zm(/) -u^v) [ NFLifu, F) = ] , (14) 

then F is called partially involutive up to the monomial v with respect to the admissible 
ordering F is still said to be partially involutive up to f if f -< lm{f) for all f & F. 

Looking at the proofs of Theorems |4.10| and |6.5| it is easy to see that they prove also 
the following conditions of partial involutivity. 

Corollary 6.7 Given a continuous involutive division L, an L— autoreduced set F is 
partially involutive up to the monomial v if and only if 

(V/ G F) (Vx, G iVMi(/m(/), lm{F))) {lm{f) ■ -< v) [ NF^if ■ x,, F) = ] . (15) 



7 Grobner Bases and Involutive Bases 



In |I3| , Q it was shown that a Pommaret basis, that is, involutive basis for Pommaret 
division, is also a Grobner basis, though, generally, not the reduced one. A similar 



property of a Janet basis was noticed in |]T^. The following theorem shows that such 
a relation holds for any involutive division. 



Theorem 7.1 // set F is L— involutive, then the equality of the conventional and 
L— normal forms 

(Vp G R) [ NF{p, F) = NFl{p, F) ] (16) 
holds for any normal form algorithm. 



Proof To prove the theorem it is sufficient to show that any polynomial p is reducible 
modulo F if and only if it is involutively reducible. But the latter statement is an easy 
consequence of Definitions |3.1| or |3]^ and |6.2| . Indeed, if p is involutively reducible, 
then it is conventionally reducible. Conversely, let the term u have a divisor among 
the leading monomials of F, that is, u = lc{u) lm{f) ■ v for some f & F and f G M. 
By the condition (^ and Theorem |5.4| , it implies f ■ v = fi ^ ^ij- Hence, u has 
also the involutive divisor in lm{F). It is just that fi which satisfies the condition 
lm{fi) X Uij = lm{f) ■ V and is unique. □ 
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Corollary 7.2 An involutive basis is a Grobner basis. 



Proof According to the algorithmic characterization of Grobner bases [|l], ^ Q consider 
the S-polynomial of fi, fj e F 

qfr ^^ _ lcm{fi,fj) ^ Icmjfijj) ^ 

itu.) im ^^^^ 



From S{fi, fj) G Id{F), Corollary U and Theorem we have NF{S{f^, fj), F) = 0. 

□ 

Corollary 7.3 // set F is partially involutive up to the monomial v, then 

(yp e R) {lm{p) -< v) [ NF{p, F) = NFl{p, F) ] . (18) 



□ 



Proof It follows by perfect analogy to the proof of Theorem 1771 . 

Note that while a Pommaret basis, if it exists for the given ideal, is unique |]14|, this 
may not hold for other involutive divisions. We demonstrate it by the following explicit 
example. 

Example 7.4 Two lexicographical {x >- y) Janet bases Fi and F2 



f \ f \ 

Fi = {xy^ -y,xy^ -\,xy -y^,x -y^y^ -1} , 

x,y 

..... . ... 



y X x,y X 

r 2 3 2 2 2 2^ i 2 2 ' s"" ' ' 2 i 2 ^ 3 7l 

{x y -y,xy -y,xy-l,x -y,xy -y,xy - l,xy - y ,x - y,y - 1} 



with indicated non-multiplicative variables, are involutive. It can easily be verified. 
Both of them generate, obviously, the same ideal with the Grobner basis {x — y,y^ — l), 
which is also a Janet basis and, in this particular case, coincides with the Pommaret 
basis. 

As it was shown in Sect. 4, given a polynomial set F and an arbitrary involutive division, 
the ideal Id{F) may not have a finite involutive basis. For example, while a finite 
Pommaret basis exists for any zero- dimensional ideal |TI|, IT3, [TH] , it may not exist for 



a positive dimensional one. Generally, for positive dimensional ideals, the existence of 
finite Pommaret basis can be achieved by means of an appropriate linear transformation 



of variables 11, 18 



On the other hand, a noetherian involutive division, for example, a Thomas or 
Janet one, implies the existence of finite involutive bases for any polynomial ideals as 
the following proposition shows. 
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Proposition 7.5 // involutive division L is noetherian, then any polynomial ideal 
Id{F) has a finite L— involutive basis. 

Proof Let G be the reduced Grobner basis of Id{F) which is finite for any polynomial 
ideal 0, If set G is not involutive, then complete it by non-multiplicative prolonga- 
tions of its elements just as it done in algorithm InvolutiveCompIetion. This means 
that at every step of the completion we select a non-multiplicative prolongation with 
the lowest leading term which is L— irreducible modulo the current leading monomial 
set. By noetherity of L, in a finite number of steps, a polynomial set G will be pro- 
duced such that lm{G) be an L— autoreduced involutive completion of lm{G). Finally, 
L— autoreduction of the tales in G will give an L-involutive basis of Id{F). □ 



8 Basic Algorithm 



In this section we describe an algorithm for the construction of an involutive basis. The 



algorithm is an improved version of one presented in for Pommaret division and 
generalized to any continuous noetherian division L and any admissible ordering y. 
The main optimization is based on the use of Buchberger's chain criterion for avoiding 
unnecessary reductions introduced in |TB[ (see also 0, 

Corollary p]3| shows that for any S-polynomial S{fi, fj), given by formula (^7^, both 
its conventional and L— normal forms are vanishing as soon as the conditions (|15]) are 
satisfied up to the monomial lcm{fi, fj). According to Theorem p.4| and Corollary p.5| 
the conditions ([TsD can be presented as NF^^SLifi, fj), F) = 0, where Si^fi^fj) are 
just (L— involutive) S-polynomials of the special form 

SlUi^ fj) = fi-x- fj X Ujk . (19) 

The following theorem gives the involutive form of Buchberger's chain criterion. 

Theorem 8.1 Let F be a finite L— autoreduced polynomial set, and let g ■ x be a non- 
multiplicative prolongation of g E F. Then NFi^g ■ x,F) = if the following holds 



(yh e F) (yu e m 
(3/,/o,^?oeF) 



{ lm{h) ■ u ^ lm{g ■ x) ) [ iVFi(/i ■ m, F) = ] , (20) 
lm{fo)\lm{f) , lm{go)\lm{g) 

lm{f)\Llm{g ■ x) , /cm(/o, 5-0) ^ lm{g ■ x) . ( 21) 
NFM-llj^^,F)=NF,igo-§^,,F)=0 



U{go) ■ 



Proof Condition (|2l|) yields that at least one of polynomials /, g can be considered as 
derived from fo,go by prolongations with at least one non-multiplicative among them. 
If, for example, /m(/o) 7^ lm{f), it leads to the equality / = /o-(/m(/)//m(/o)) modulo 
F. 
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Thus, if the condition ([21| ) holds, there is a chain of polynomials in F of the form 



/ = /fc) fk-l-, ■ ■ ■ , fo, 90, ■ ■ ■ , Qm-l, dm = 9 , (22) 

where k + m > 0. Here f or g or both of them are produced by prolongations, including 
non-multiplicative ones, of the polynomials fi or gj in the chain whose indices are less 
than k or m, respectively. 

The chain (|22|) has the property 

NF{SL{f, fk-i), F) = ■ ■ ■ = NF{S{fo, go),F) = --- = NF{SL{g„.-i, g), F) = . 

This property is resulted from the observations as follow. Consider relation 

lm{g) ■ X = lm{f) x w, (23) 

which means that w does not contain x. Otherwise, g would be reducible by /, and, 
hence, F could not be L— autoreduced. Thus, lcm{f,g) = lm{g) ■ x. By admissibility 
of the monomial ordering -<, the least common multiple of the leading monomials 



for pair of the neighboring polynomials in the chain (22) is less than or equal to g ■ x. 
Then the above property of the chain follows immediately from partial involutivity (|20|) 
of F and Corollary [f73[ Furthermore, conditions (pO| - pTD imply NFL{S{fo, go), F) = 
NFiS{fo,go),F) = 0, and NFL{SL{f^, f^-l), F) = NF{S{f,, f,^i), F) = (1 < t < k) 
as well as NFL{SL{gi-u gi), F) = NF{S{gi.i, gi), F) = 0{l<t<m). 

By construction, lcm{f, . . . , /i, /q, go, gi, . . . ,g) = lcm{f, g) what leads to the 
representation S{f,g) = fiUij where fi & F and lm{fiUij) -< lcm{f,g) = lm{g) ■ x. 
Then, condition (^0]) , by Corollaries and [773| , yields 

NFl^SlU, g),F) = NF{S{f, g),F) = 



in accordance with [Q, 0. □ 



Before analysis of correctness and termination of the below algorithm, we give some 
necessary clarifications. 

First of all, the conventional autoreduction of the initial polynomial set is done. It 
removes, in particular, all the predecessors of every polynomial from the initial set. 

Set T collects all the triples {g,u,P); g is an element in the current basis G; u = 
lm{f) where / G G is the predecessor of g, by a non-multiplicative prolongation of 
which g was derived, or m = lm{g) if g has no such predecessor in G; P is a set 
containing the non-multiplicative variables of g have been used for its prolongations. 

The current non-multiplicative prolongation • x is selected to be the lowest with 
respect to the ordering >-. If there are several different non- multiplicative prolongations 
with the same leading term, then any of them may be selected. This selection strategy 
will be called normal. 

If the leading monomial of the current prolongation g-x is involutively reducible by 
the basis element f & G, then the other conditions in (plf) are verified. The verification 
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is done in the form of comparison of lcm{u,v) with lcm{f,g), where u and v are the 



second elements of the triples containing g and /, respectively. By Theorem the 
criterion (^Tj) is false if and only if lcm{u,v) = lcm{f,g) = g ■ x. One should be also 
noted that Buchberger's second criterion can be applied in the involutive approach 
only in exceptional cases. Relation ( p3D shows that Zcm(/, (7) = lm{f)lm{g) if and only 
if lm{f) = X and lm{g) = w. 

If the current prolongation is not reducible to zero, that is, h = NFi^g ■ x,G) 7^ , 
then h is added to G. 

After involutive autoreduction of the enlarged set G an adjustment of the set T 
is done. For an element g & G whose leading monomials was not mutually reduced, 
the second element u in the triple is kept, if the leading term of the corresponding 
predecessor of g was also not reduced. Otherwise, u is replaced by its involutive divisor 
in lm{G). Essentially new leading monomials, that is, those not multiple of any others 
occurring in T before the autoreduction, are included in the refreshed T with their 
actual leading monomials as the second elements of the triples. 

Algorithm InvolutiveBasis: 

Input: F, a finite polynomial set 

Output: G, an involutive basis of the ideal Id{F) 

begin 

G := Autoreduce{F) 
T := 

for each g e G do T := T U {(^, /m(^), 0)} 

while exist {g, u,P) eT and x G NML{lm{g), lm{G)) do 

choose such {g, u, P) and x with the lowest lm{g) ■ x 

T:=T\{{g,u,P)}VJ{{g,u,PVJ{x})} 

if exist / G (/, v,D) & T such that lm{f ) \l lm{g ■ x) then 
if lcm{u,v) = lm{g) ■ X then h := NFilg ■ x,G) 
ii h^O then T := T U {(/;., /m(/i), 0)} 
else h := NFL{g ■ x,G) 
T:=TU{(/i,u,0)} 
G := AutoreduceL{G U {h}) 
Q:=T 
T := 

for each g E G do 

if exist {f,u,P) G Q such that lm{f) = lm{g) then 

choose gi E G such that lm{gi)\LU 

T:=TU{{gM{9i).P)} 
else T ■=TU{{gM{9)M 

end 
end 

Correctness. As we have shown, criterion ( ^1]) is used in algorithm InvolutiveBasis 
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in accordance with Theorem It is easy to show that there is the unique polynomial 
gi ^ G which is chosen in the inner for each-loop such that lm{gi) involutively divides 
u. Indeed, if the leading term of the predecessor h of g with u = lm{h) has not been 
reduced, then gi = h. Otherwise, there is gi E G such that gi ^ h and lm{gi)\LU. The 
uniqueness of gi for the autoreduced set G is an immediate consequence of the property 
(v) in Definition |3.1| . Besides, the replacement of u by gi does not violate, obviously, the 
conditions for applicability of the criterion. Furthermore, from Corollary |7.3| it follows: 
(i) a leading monomial, being involutively reducible at some step of the algorithm, will 
never appear again among the leading monomials; (ii) there is no need in recomputing 
zero reductions after enlargement of an intermediate polynomial set. This enables 
one to assign the set P of the used non-multiplicative variables for polynomial / to 
the corresponding polynomial g with lm{g) = lm{f) as it is done in the inner for 
each-loop. Such an optimization allows one to avoid the repeated prolongations. 

Therefore, if division L is continuous, and the algorithm terminates, then it pro- 
duces, by Theorem |6.5| , the involutive basis. The termination holds if and only if the 
set P in each triple {g, u,P) E T contains all non-multiplicative variables for basis 
element g. It just means that any non-multiplicative prolongation of every element in 
G is reduced to zero, and, hence, G is involutive. 

Termination. Note that the initial value of the leading monomial set 

Uq = lm{Autoreduce{F)) 

is determined by the input set F subjected to the conventional autoreduction. Since 
only those monomials occur in the leading monomial set which have not been reducible 
at some step of the algorithm, the change in set U = lm{G) after running the while- 
loop may take place only in two cases: 

(i) . lm{g) ■ X has no involutive divisors in U . In this case U is enlarged to include 

lm{g) ■ X. 

(ii) . g ■ X is reducible by elements of U. Then U is enlarged to include lm{h), where 

h = NFi^g ■ x,G) 7^ and lm{h) is not multiple, in the conventional sense, of 
any elements in f/o. 

The number of different lm{h) occurring in case (ii) is finite by Dickson's lemma 
(Becker, Weispfenning and Kredel, 1993). Recall also that algorithms InvolutiveAu- 
toreduction and InvolutiveNormalForm always terminate (Sect. 5). 

Thus, the algorithm termination is determined by that of algorithm Involutive- 
Completion considered in Sect. 4. It follows that algorithm InvoIutiveBasis termi- 
nates for any noetherian division and arbitrary input polynomial set F . If L is not 
noetherian, then termination may not hold if an intermediate set U = lm{G) is not 
finitely generated with respect to L as the below Example ^.2| shows. In the case 
of Pommaret division the algorithm terminates, however, for any degree compatible 



ordering and any zero-dimensional ideal WM. Because the involutive division L is con 



tinuous, once algorithm InvolutiveCompIetion terminates, an L— completion U of 
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U will be constructed such that autoreduction of the corresponding set G does not 
produce new leading monomials. G is, obviously, the output involutive basis. 

Proposition implies, in particular, the algorithm termination for Thomas and 
Janet divisions. However, for Pommaret division, which is not noetherian, the algo- 
rithm may not terminate even in the case when there is a finite Pommaret basis but 
the ordering is not degree compatible as the following simple example shows. 

Example 8.2 The set F = {x^ — l,xy — l,z} generates a zero-dimensional ideal with 
the lexicographical Pommaret basis {x >~ y >~ z) given by G = {x — y,y'^ — l,yz,z}. 
However, following the above algorithm we have to choose z-y a.s the first prolongation 
which is lexicographically lowest. Since polynomial h = yz has no Pommaret divisors 
among lm{F), we find FU{yz} as an intermediate basis. The next lowest prolongation 
is yz -y again has no Pommaret divisors among the leading monomials of the enlarged 
set. Exploring this procedure further produces the infinite involutively irreducible set 

{x^ -l,xy- 1, z, yz, y^z, . . . , y'^z, . . .} k eN. 



It is well-known [L3| , ^4| , [TSf that positive dimensional ideals may not have finite 
Pommaret bases. Example |4.16| illustrates this fact at the monomial level. The fol- 
lowing more non-trivial example shows the output of algorithm InvolutiveBasis for 
Pommaret and Janet divisions in the case of polynomial ideal. 

Example 8.3 Cyclic 4-th roots. 



NMj 


NMp 


Initial Polynomial Set 






Xi + X2 + Xs + Xi 




Xi 


X1X2 + X2XS + X3X4 + X4X1 


X4 


Xi,X2 


X1X2XS + X2X2,Xi + X3X4X1 + Xi^XiX2 




Xi,X2, X3 


X1X2X3X4 — 1 



Here we choose the degree-reverse-lexicographical-ordering with the order of variables 
as in Sect. 2. Note that, since the initial set is not autoreduced, the inclusion NMj C 
NMp (see Proposition p.lOj ) does not hold. 

Application of algorithm InvolutiveBasis gives the following form of Janet and 
Pommaret bases 
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NMj 


NMp 


Janet and Pommaret Bases 












Xi + X2 + X3 + X4 




Xi 






Xi 






Xi,X2 




Xi 


X2 


X2xl + X3X4 — X2xl — xl 


Xi 






To 


To 


a/2'^3'^4 '^3'^4 ■^2'^4 ~r ^S^^^ "^4 


Xi 




Xi 


X2 


X3 


X2x\ + x| — X2 — X4 


Xi 


X2,X3 


Xi 


X2 




X3X4 + X2X3 — X2X4 + X3X4 — 




Xi,X2 


Xi 


X2 


X3 


3 4 ~'~ '^3*^4 "^4 




Xi 


X2 


X3 


x|x4 + X2X3 — 0:3 — X2X4 + X3X4 — X4 



The Janet basis consists of the upper seven polynomials and coincides with the Grobner 
basis, while the Pommaret basis is infinite and contains also prolongations of the sev- 
enth polynomial with respect to its non- multiplicative variable 0:3. Note that the ideal 
is one-dimensional, that is why it does not have a finite Pommaret basis. 

The algorithm InvoIutiveBasis has been implemented in Reduce 3.5 for the degree- 
reverse-lexicographical-ordering and Pommaret division refined in a certain way to pro- 
vide the algorithm termination for any polynomial ideal. This refinement is equivalent 
to the dynamical incorporation of some noetherian involutive division in the compu- 
tational process. Its detailed description will be given elsewhere. In addition, the 



current package called INVBASE is considerably faster than previous version |jTj], in 
particular, since it uses the criterion (pTf). 

Experimentally, we observed much smoother behavior of the algorithm Involutive 
Basis with respect to Buchberger algorithm^ as the ordering changes. Consider, for 
instance, the following example. 



Example 8.4 Cyclic 6-th roots. 

Xl X2 + X3 + X4 + X5 + , 

X1X2 + X2X3 + X3a;4 + 0:4X5 + X5X6 + xeXi 



X1X2 + X2X3 + X3X4 + X4X5 + X5X6 + xeXi , 
X1X2X3 + X2X3X4 + X3X4X5 + X4X5X6 + X5X6X1 -f X6X1X2 , 
X1X2X3X4 + X2X3X4X5 + 0:3X4X5X6 + X4X5X6X1 + X5X6X1X2 + X6X1X2X3 , 
X1X2X3X4X5 + X2X3X4X5X6 + X3X4X5X6X1 + X4X5X6X1X2 + X5X6X1X2X3 + X6X1X2X3X4 , 
X1X2X3X4X5X6 — 1 . 

The next table gives the timings of INVBASE on an 66 Mhz MS-DOS based AT/486 
computer for different degree-reverse-lexicographical-orderings. 



^More precisely, with respect to its implementation in Reduce 3.5. 
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Ordering 


Timing (sec.) 


Xi 


y 


X2 


y X3 


y 


X4 


y 


X5 


y 


Xe 


1040 


Xi 


y 


X2 


>- X4 


y 


Xe 


y 


X3 


y 


X5 


514 


Xi 


y 


X2 


>- X4 


y 


Xe 


y 


X5 


y 


X3 


437 


Xi 




X2 


y Xe 


y 


X3 


y 


X4 


y 


X5 


1066 


Xi 






>- Xi 


y 


X5 


y 


X2 


y 


Xe 


604 


Xi 


y 


X3 


>- Xi 


y 


Xe 


y 


X5 


y 


X2 


136 


Xi 


y 


X4 


y X2 


y 


X3 


y 


X5 


y 


Xe 


993 


Xi 


y 


X4 


y X5 


y 


Xe 


y 


X2 


y 


X3 


1001 


Xi 




X5 


y X3 


y 


X4 


y 


Xe 


y 


X2 


364 


Xi 


y 


X5 


y Xe 


y 


X2 


y 


X3 


y 


X4 


1045 


Xi 


y 


xe 


y X3 


y 


X2 


y 




y 


X5 


1012 


Xi 


y 


xe 


y X5 


y 


X2 


y 


X4 


y 


X3 


590 



Comparison with the package GROEBNER implementing Buchberger algorithm on 
the same Reduce 3.5 platform shows that its corresponding timings are not only much 
larger than those presented in the table, but also vary dramatically with the order of 
the variables. This fact was already observed in |jT^ where some comparative data for 
GROEBNER and the previous version of the INVBASE package are presented. 



9 Conclusion 

Buchberger algorithm and the involutive one are based on different rewriting tech- 
niques, namely, on the use of S-polynomials and prolongations, respectively, as well as 
on distinct reduction processes. Nevertheless, as we demonstrate in this paper, they 
are in fact very interconnected. If, as we propose in the algorithm InvoIutiveBasis, 
we choose the current prolongation in increasing order with respect to given monomial 
ordering, then the conventional and involutive normal form will coincide. What is 
more, the involutive reduction of the prolongation is equivalent to the consideration of 
a certain S-polynomial. Just this fact makes it possible to use Buchberger's criteria. 

Recently another interesting facet of interrelation of both methods was discovered 
by Apel |T^, namely, that Pommaret bases can be associated with Grobner ones in 
appropriate graded structures. Earlier such Grobner bases were intensively investi- 
gated in more general context by Mora [Q. That observation gives an opportunity to 
algorithmically construct Pommaret bases whenever they exist |1^. Though such an 



analogy also enables one to take advantage of Buchberger's criteria, it is restricted to 
Pommaret division. 

Thus, all the above, as well as computer experiments with both techniques, offers 
a clearer view of the most optimal computational procedures. 

There is no question that any algorithmic improvement of the Grobner basis and 
involutive techniques at the algebraic level has an analogous optimization at the dif- 
ferential level, at least for linear partial differential equations [0 . 
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